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I. INTRODUCTION 

The quantum Heisenberg model (QHM)0 and its simplified Ising versioni remain to be one of the most 
actively studied subjects in statistical mechanics. In particular, the low-dimensional antiferromagnetic 
QHM has recently attracted a lot of attention since it represents a very useful model for the investigation 
of interesting quantum phenomena. Among the most fascinating problems that are of current interest in 
this field oncj-should mention: the investigation of Haldane gaps in the. ID antiferromagnetic QHM with 
integer spiiisu, quantum phase transijtionsu, spin-Peierls instabilitiesEI, dimerization and other related 
phenomenal^, quantum entanglement)]], magnetization plateauu, and so on. In addition to the above 
mentioned works, a number of the studies has been devoted to the investigation of the role of magnetic 
ordering in high-T c superconducting cuprates consisting of two-dimensional Cu-0 networksa. In fact, the 
magnetic properties of these materials can be well-described by means of a spin-1/2 antiferromagnetic 
2D QHM. _ 

Despite of extensive studies, the QHM has been exactly solved in one dimension onlyE3, thus for 
higher dimensions only more or less accurate approximate methods are available^. In general, the main 
difficulties of a rigorous treatment of the QHM are closely associated with the noncommutability of the 
spin operators involving the Hamiltonian of the system. This principal mathematical intractability of 
the QHM has motivated us to introduce a class of interesting exactly solvable models, consisting both of 
the Ising- and Heisenberg-type atoms. For this purpose, we utilize the well known decoration-itecation 
procedure that has been originally introduced by SyoziO and later remarkably generalized by Fisheniij. In 
the spirit of the Syozi's and Fishers's papers, we use in this work the decoration procedure to put a couple 
of the Heisenberg atoms on each bond of the regular Ising lattice. In this way we obtain a doubly decorated 
model consisting of two interpenetrating sublattices occupied by the Ising- and Heisenberg-type atoms, 
respectively, and such a model can be of interest both theoretically and experimentally. Theoretically, the 
most outstanding feature of this model is the fact that it enables to investigate, by an exact calculation, 
how the quantum Heisenberg atoms modify the magnetic properties of the pure Ising systems. On the 
other hand, from the experimental point of view, the model can be inspiring in the preparation of new 
magnetic materials with similar topological structure as the system under investigation. Actually, some 
of recently synthesized compounds represent a progressive step in this directions. 

The outline of the present paper is as follows. In Sec. II, the main points of the mathematical 
formulation of the decoration-iteration procedure for the Ising-Heisenberg models are explained and the 
exact equation for the phase diagrams and basic physical quantities are derived. The most interesting 
numerical results are presented and discussed in detail in Sec. Ill and finally, some concluding remarks 
are given in Sec. IV. 
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II. FORMULATION 



In this work we will study a spin- 1/2 doubly decorated Ising-Heisenberg model on planar lattices 
described by the Hamiltonian: 

H d = - £ J[A(S?S? + SfS*) + S?S?] - JiHfi, (1) 

i,j k,l 

where Sf and (if are the well-known components of spin-1/2 operators, the parameters J and J\ rep- 
resent the nearest-neighbor exchange interactions that couple the pairs of Heisenberg atoms or Ising 
and Heisenberg atoms, respectively (see Fig. 1). As usual, A describes the spatial anisotropy in the 
Heisenberg exchange interaction. In fact, this parameter allows one to control the behavior of the system 
between the Ising-regime for A < 1 (easy-axis like anisotropy) and the XY-regime for A > 1 (easy-plane 
anisotropy). In view of further manipulations, it is useful to rewrite the total Hamiltonian of the system 
in the form Hd = ~^2 k Hk, where the summation is carried out over all the pairs of Heisenberg atoms. 
The bond Hamiltonian Hk involves all the interaction terms associated with fcth couple of Heisenberg 
atoms (see Fig. 1) and it is given by 

H k = -j[A(S* kl S x k2 + S y kl S y k2 ) + S* kl S* k2 ] - Ji(&£jEi + S| 2 £I 2 ). (2) 

The most important point of our treatment is the calculation of the partition function for the system 
under investigation. Taking into account the standard commutation relation for the bond Hamiltonians 
(i. e., [Hi, Hk] = 0, i^i), we can express the partition function Zd of decorated system in the form 

Nq/2 

Z d = Trexp(-/3W d ) = Trexp(-/5 ^ H k ) 

k=l 

Nq/2 

= Tr^HTrsJ^exp^), = l/k B T, (3) 
fc=i 

where ks being Boltzmann constant and T the absolute temperature. N represents the total number of 
Ising atoms and q is the coordination number of the original (undecorated) lattice. The symbol Tr{ M } 
means a trace over all degrees of freedom of Ising spins and finally, Tr5 )tl Tr5 fc2 denotes a trace over a 
couple of Heisenberg spins residing on the kth bond. _To_proceed further, it is useful to introduce the 
following extended decoration-iteration transformatior£3J13 

Tr Sfcl Tr 5fc2 exp(-/37^) = 2exp(/3J/4)|cosh[/3J 1 (^ 1 + fi z k2 )/2 

+ exp(-{3 J/2) cosh [^y/jR^x - M^ 2 ) 2 + J 2 ^ 2 ] } = Aexp((3R^ z k2 ). (4) 

The unknown transformation parameters A and R can be easily obtained following the standard proce- 
dures (seelia and references therein), namely, 



A = 2exp((3J/A)(V 1 V 2 ) 1 ^ 2 , /?i? = 21n(^), (5) 



' V ± 

where we have introduced the functions V\ and V 2 as follows 



Vi = cosh(/3Ji/2) + exp(-/3J/2)cosh(/3JA/2), 

V 2 = 1 + exp(-/?J/2) cosh(|^J? + J 2 A 2 ) . (6) 
Now, after substituting Eq. (4) into Eq. (3) one obtains the equation 

Z d = A N ^ 2 Z , (7) 

which relates the partition function of doubly decorated Ising-Heisenberg model {Z d ) to that of the 
original undecorated spin-1/2 Ising model (Zq.) From this simple relation, we can directly calculate 
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some physical quantities (for example, the free and internal energy or specific heat) on the basis of well- 
known thermodynamic relations. However, to understand the behavior of the system, we have to analyze 
also some other quantities (for instance the magnetization and pair-correlation functions) that cannot 
be obtained from the partition function in a straightforward manner. Fortunately, we can avoid this 
complication by exploiting the following exact identities 



(/i(/tf,A?,->Afe))rf = (/i(Az,A?,---,AD)o, 



/ Tr Skl Tr Sk2 f2(Ski , S% 2 , AjjjAfg) exp(-(3H k ) 
\ Tr Sfcl Tr S)c2 exp(-/37i fc ) 



(8) 
(9) 



from which the relevant quantities can be calculated. In the above equations, f± represents a function 
depending only on the Ising spin variables and fa denotes a function which depends on the spin variables 
located on the fcth bond only. The superscripts a and 7 denote x, y or z components of the spin operators 
and the symbols (...)<j and (...)o mean the standard ensemble average related to the decorated and original 
lattice, respectively. Now, applying one of the standard methodst3, we simply derive the following 
equations for sublattice magnetization 



m A = ^((Afci + Afc2>d) = ^((Afci + Afc2)o) = mo, 
m z B = ^(S^ + S z k2 ) d = 8«Afci>d + (fik2)d)K 



(10) 

(11) 
(12) 



■.y — 



S 



I2U = 



0. 



(13) 



where mo is the magnetization per one site of the original lattice and the coefficient Kq is given in the 
Appendix. Similarly, the various pair correlations can be expressed with the help of Eqs. (|^) and (|9|) in 
the following simple form: 



XX 

Ihh 



tiri 

Qhh 



Qhh 



(§kA)d = K 1+ K 2 + 4qff(Kt - K 2 ), 

Wfel^D^ = Ihhi 

(S z kl S z k2 ) d = K 3 + K i + Aq zz (K 3 - K 4 ), 



-(S^AJIi + S z k2 fi z k2 ) =K + K 5 + Aq zz (K - K 5 ). 



(14) 
(15) 
(16) 
(17) 

(18) 



Here, e = (AfciAfc2)o denotes the nearest neighbor correlation of the original lattice that is well known 
and the coefficients Kq — K$ are listed in the Appendix. 

Finally, the internal energy and specific heat of the system can be also easily calculated from the 
relations 



U d 



Nq 
~2 



jA(qr h + qyy h ) + j q zz h + 2j iq zz 



(19) 



c d = dujdT. 



(20) 
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III. NUMERICAL RESULTS AND DISCUSSION 



In this section we will show the most interesting numerical results of the system under investigation. 
For the sake of simplicity, we restrict our attention to case of the doubly decorated square lattice (see 
Fig. 1) in which all characteristic properties can be illustrated. 

Before discussing the results, it is worth noticing that the phase diagrams for the ferromagnetic (J > 
0, Ji > 0) and ferrimagnetic (J > 0, J\ < 0) case will be the same, since the relevant equation for 
the critical temperature is invariant under the transformation J\ < — > —Ji- On the other hand, the 
antifcrromagnetic system ( J < and arbitrary Ji) exhibits many different features and will be discussed 
in a separate work. 

In order to find possible ground state phases and to investigate their properties, we have to analyze the 
internal energy, magnetization and correlation functions at T = 0. Depending on whether the anisotropy 
parameter A is less, equal or greater then the boundary value A c = y/2\ J\\/ J + 1 one finds three different 
regions. Namely, 



i) for A < A 

Nq 



= 0.25, g « = ±0.25 > qH = qZ = G.O vll) 



ii) for A > A c 

U d = -^(-J + 2yJj* + (JAy 
C=Cfc = -0-25, q^ = ± ' h 



VJ! 2 + (JA)2' 

VJ? + (JA)2 



f h l = q ™ = —JJt== (22) 



iii) for A = A c 



2\J\\ + J xx _ yy _ yj J{\Jl\ + J) 



q" -4l-o.o, C-± 8(|Ji| + J); €l-% h - 8(1^1 + j) (23) 

where the plus or minus sign one applies for the ferromagnetic or ferrimagnetic case, respectively. From 
these relations, we have obtained the ground-state phase diagram in the |Ji| — A space which is de- 



picted in Fig. 2. Taking into account Eqs. (]2ll)-(23), one easily identifies the standard ferromagnetic 



(ferrimagnetic) phase (FP) for A < A c . However, for A > A c an unexpected quantum phase occurs 
in the system. This phase (to be referred to as a quantum antiferromagnetic phase (QAP)) requires 
more detailed description since it differs from standard phases in the pure Ising or Heisenberg models 
and, as far as we know, such a phase has not been described in the literature before. In fact, from 
the relevant equations we find that in the QAP the Ising spins (that are the nearest neighbors on the 
original lattice) are aligned antiparallel with respect to each other. Consequently, we have the classical 
Neel long-range ordering on the Ising sublattice with uia — 0.0 and qff = —0.25. On the other hand, 
the nearest-neighboring Heisenberg spins create dimers, thus we also have tub — 0.0 and qfX = —0.25. 
However, the alignment of these dimers with respect to their nearest-neighboring Ising spins is random, 
hence, the relevant correlation function qf? does not reach its saturated value (±0.25). Moreover, one 
easily observes that the degree of randomness increases with the increasing in anisotropy A. This be- 
havior apparently appears due to the competition between strong in-plane anisotropy (A) that supports 
short-range ordering and the exchange interactions (J, J\) preferring the long-range ordering along the 
easy axis. It is also clear from the aforementioned arguments that despite of some disorder introduced 
by random orientation of the dimers, the QAP will exhibit the long-range Neel ordering captured to the 
Ising spins. Thus, one can expect the appearance of the second-order phase transition in the system at 
finite temperatures, even for very strong values of the anisotropy parameter A. Nevertheless, the QAP 
also differs from the standard Neel phase (both the classical and quantum one) due to the ferromagnetic 
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in-plane short-range order of the Heisenberg spins. Another unusual feature of the QAP is the perfect 
antiparallel alignment of the Ising spins (that are not directly coupled via exchange), despite of some 
disorder present between the nearest-neighboring Ising and Heisenberg atoms. Furthermore, we would 
like to emphasize that the existence of the QAP by itself is very surprising, since we have the system 
with the ferromagnetic exchange interactions J and Ji, only. 

Finally, one should notice that the FP and QAP are se parated by the first-order phase transition line 
that is given by the condition A = A c = y/2\ J\\/ J + 1. At arbitrary point of this line there coexist 
two of the above mentioned phases (FP, QAP) together with a disordered phase (DP) in which we have 
m z A = m z B — qff — qfi — 0.0 and the nonvanishing short-range ordering both in the xy plane and along 
the easy axis (q%% — qf^ ^ 0, q?% ^ 0). The coexistence of these three phases follows from the fact that 
the relevant ground-state energies take the same value (in fact, lim A ^ A - Ud(FP) = lim A _ >A + [/^(QAP) = 

J7d(DP) A =Aj. 

Next, in order to demonstrate the overall dependences of the correlation functions on the anisotropy 
parameter A, we have depicted in Fig. 3 the relevant pair correlations for the case Ji/J = 1.0. In 
agreement with the arguments given above, we find the FP for A < and the QAP for A > V3. 
Moreover, Fig. 3 indicates that for A > \^3, the correlations q?f decreases with the increasing A and 
vanishes in the limit of A — ► oo. On the other hand, the correlation qff increases with the parameter A 
and approaches its saturation value for A — > oo. Contrary to this behavior, the correlations qff and qii 
take the saturation value independently of A, excepting the special point A = A c where they jump to 
zero. Thus, the antiparallel orientation of the relevant Ising and Heisenberg spin pairs is not affected at 
all even by the very strong anisotropy. Finally, to complete the ground-state analysis, we have shown in 
Fig. 4 the internal energy of different phases as a function of the anisotropy parameter A for | J\\/J = 1.0. 
In this figure, the full and dashed lines represent the stable and unstable parts of the relevant energies, 
respectively. The black point is the energy of the DP. This dependence apparently supports our previous 
statements and clearly illustrates the occurrence of the first-order phase transition at A = A c = \/3. 

Now, let us proceed to study the finite-temperature phase diagrams for q — 4 that can be easily 
obtained by putting j3 c R = ±21n(l + \/2) into Eq. (|^). Solving the relevant equation numerically for 
some characteristic values of J%, we have obtained the phase boundaries in the A — T c plane that are 
plotted in Fig. 5. As one can see, the critical temperature decreases gradually from its Ising value 
at A = and vanishes for A = A c = y / 2\Ji\/J + 1. On the other hand, for A > A c the transition 
temperature at first rapidly increases, then passes through a local maximum value and finally tends to 
zero for A — > oo. It is clear that for A < A c the phase boundary separates the FP and the DP, and 
similarly for A > A c separates QAP and DP. As one can expect, the relevant thermal phase transition 
is of the second order and belongs to the same universality class as that of the usual 2D Ising model. Of 
course, the thermal variations of physical quantities can differ from the standard behaviors in the Ising 
model. To illustrate the case, we have shown in Fig. 6 the temperature dependences of the specific heat 
for Ji/J = 0.5. As we can see, in the isotropic case (A = 1.0) we have the standard dependence usually 
observed in the Ising models. On the other hand, for the values of the exchange anisotropy close to 
the critical value (A c = \/2), the specific heat may in addition to the familiar Schottky-type maximum 
exhibit another maxima. These maxima appear equally bellow (see the case A = 1.3) and above (see the 
insert in Fig. 6) the critical temperature, as a consequence of the thermal excitations of the Heisenberg 
spins that basically influence the ordering in the system. The influence of the thermal excitations on 
the behavior of the system is really of great importance and can be understood from the temperature 
dependences of the correlation functions. For this purpose, we have plotted in Fig. 7, 8 and 9 thermal 
variations of the correlations qffi, q^ z h and q%? for the same values of J\ and A as in Fig. 6. It is clear from 
these figures (see the case A = 1.3) that the occurrence of the maximum bellow the critical temperature 
require to satisfy two conditions: i) the relevant correlations must increase (decrease) rapidly enough 
with the increasing temperature in the low-temperature region (for example, in the case of A = 1.0 the 
additional maximum does not appear due to the relatively slow excitation process in comparison with 
the case of A = 1.3). ii) the critical temperature of the system must be relatively high (for instance, in 
the case of A = 1.4 we have &bT c /J ps 0.0075 that is very low for the occurrence of the maximum below 
T c , in spite of the very strong thermal excitations). 

In addition to this behavior, the double-peak specific heat curve can be also observed here. The origin 
of the relevant maxima above the critical temperature can be understood from the thermal dependences of 
the correlation functions, as it is clearly displayed in Fig. 6-9 for A = 1.42. As one can see, the correlation 
function qf? rapidly increases in the relevant region, though the short-range ordering of Heisenberg atoms 



5 



in the xy plane as well as in the z direction is rapidly destroyed. Very similar behavior appears above T c 
also for A < A c excepting the fact that in this case the correlations in the xy plane abruptly increases, 
although all the other correlations rapidly decreases as it is apparent from Figs. 7-9 for the case of 
A = 1.4. Hence, the appearance of the multiple peaks in the specific heat curve arises due to the relevant 
thermally induced short-range ordering or disordering in the system. Finally, one should also mention 
that all pair-correlation function exhibit at critical temperature weak energy-type singularity known from 
the usual Ising models. 

Although, the numerical calculation have been presented for the square lattice q = 4, we can on 
the basis of our formulation draw some general conclusions about the behavior of the doubly decorated 
Ising-Heisenberg systems. However, one should emphasize that our next statements do not concern the 
one-dimensional case that exhibits quite different behavior, as result of the fact that no long-range order 
is possible in the system at nonzero temperatures. 

Firstly, it is clear from Eqs. (f2lj)-(p3"|) that the ground-state phase diagram does not depend on the 
coordination number and spatial dimensionality of the system. This implies that the value of A c is also 
independent of the coordination number and dimension, although our preliminary investigation of the 
other systems has revealed that it depends on the spin value of the Heisenberg atoms. In general one can 
say that the QAP exists in many two- and three-dimensional lattices in the region of A > A c . However, 
it is necessary to emphasize the fact, that the above statement is valid only in the case when the Neel 
order is possible on the relevant original lattice. Apparently, this not the case, for example, of the doubly 
decorated triangular or Kagome lattices for which much more complicated phases will occur for A > A c . 

It is also interesting to note that the QAP, as well as the standard FP may exist at finite temperatures 
and the temperature region of their stability is clearly enlarged with the increasing coordination number 
and dimensionality of the system. This is a consequence of the fact that the critical temperature increases 
with the increasing coordination number and dimensionality of the system. Of course, the thermal 
fluctuations gradually destroy the long-range order (both for the FP and QAP) in the system and if 
the temperature reaches its critical value, then the system undergoes the second-order phase transition 
sharing the same universality class as the standard spin- 1/2 Ising planar model. 

IV. CONCLUSION 

In this work we have studied the doubly decorated Ising-Heisenberg model on planar lattices. Ap- 
plying the standard decoration-iteration transformation, we have obtained the simple relation between 
the partition function of the decorated model and its corresponding standard spin-1/2 Ising model. On 
the basis of this relation, we have derived the exact results for basic physical quantities that have been 
subsequently discussed for J > in Sec. III. 

In our opinion, the prediction of the QAP is the most important finding of this work and it illustrates 
how the magnetic properties of the pure Ising systems can be essentially modified by introducing the 
quantum Heisenberg atoms on the bonds of the original lattice. The origin of the QAP consists in the 
quantum fluctuations arising in the system, hence the QAP itself provides a clear manifestation of the 
quantum phenomena in the macroscopic scale. 

Although at the present time we are not aware of any experimental system which can be directly 
described by the considered doubly decorated Ising-Heisenberg model, we hope that the recent progress 
in molecular engineering will result in the preparation of such materials. In fact, recently synthesized 
compound K3C0 (CN)g.2Rhp(CH3COO)4 (see Ref. [16]), which has the magnetic structure of the 
doubly decorated square lattice (Fig. 1), seems to be the most promising from this point of view. 
Unfortunately, the Co 3+ ions (located at the corners of each square) are in this compound due to the very 
strong ligand field of the cyanide groups in the low-spin state (i.e. they are diamagnetic) . Nevertheless, 
tuning the ligand field around the Co 3+ ions by the choice of other ligands, represents a possible way 
how to prepare the compounds with high-spin Co 3+ paramagnetic ions. Another possibility to obtain the 
compounds of desired magnetic structure consists in the chemical replacement of the Co 3+ ions by other 
transition metal ions, such as Fe 3+ and Cr 3+ ions. Because the corner atoms of each square would possess 
spin 1/2 (in case of Fe 3+ ion) or 3/2 (for Cr 3+ ion), the present theory could be applied to describe the 
behavior of these materials. Since the most of the real materials have A w 1, it is also worth noticing 
that in the case of very weak exchange interaction J± (| J±\ << J), the QAP appears in our system near 
A s» 1 regardless of the coordination number and spatial dimensionality. This supports our hope that 
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the experimental confirmation of this phase would be possible in some materials. 

Finally, we would like to remark that the present formalism can be extended to investigate many other 
interesting systems. Indeed, we have succeeded in solving some interesting generalization of the system 
studied in this paper and we have found very rich and interesting behavior that will be discussed in 
forthcoming works. 

Acknowledgment: This work has been supported by the Ministry of Education of Slovak Republic 
under VEGA grant No. 1/9034/02. 



Appendix 

1 sinh(/3Ji/2) 



K = 
K x = 

K 2 = 

K 3 = 

Ki = 

K* = 



8 cosh(/3J!/2) + exp(-/3J/2) cosh(/3JA/2) ' 

1 sinh(/3JA/2) 

8 cosh(/3JA/2) + cxp(/3J/2) cosh(/3Ji/2) ' 



1 J A sinh(^/J 2 + J 2 A 2 /2 

8 y/J? + J 2 A 2 cosh ^^j2 + j2 A 2/2) + cxp(/3J/2) ' 

1 cosh(/3Ji/2) - exp(-/3J/2) cosh(/3JA/2) 
8 cosh(/3Ji/2) + exp(-/3J/2) cosh(/3JA/2) ' 

1 cxp(/3J/2) - cosh(/3 v /j 1 2 + J 2 A 2 /2) 

8 exp(/3J/2) + cosh(/3y / J 1 2 + J 2 A 2 /2) ' 

1 j, sinh(/?v^? + J 2 A 2 /2) 

8 v / jfT7 2 A 2 cosh(/3 v / JFTJ 2 A 2 /2) + exp(/3J/2) ' 



The detailed derivation of the coefficients given above requires lengthy calculation, thus is not presented 
here but it can be obtained by the authors on the request. 
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Figure captions 



Fig. 1. The fragment of the doubly decorated square lattice. The black (grey) circles denote the positions 
of the Ising (Heisenberg) atoms, respectively. The ellipse demarcates a typical bond described by 
the Hamiltonian Ttk (see Eq. (2)). 

Fig. 2. Ground-state phase diagram of the doubly decorated Ising-Heisenberg model. The full lines rep- 
resents the line of first-order phase transitions that separates the ferromagnetic or ferrimagnetic 
phase (FP) from the quantum antiferromagnctic phase (QAP) 

Fig. 3. Dependences of the pair-correlation functions at the frond state on the anisotropy parameter A for 
the doubly decorated square lattice (q = 4) and | Ji\/J = 1.0. 

Fig. 4. Ground-state energy Ud v.s. A for q = 4 and \Ji\/J = 1.0. The full and dashed lines represent, 
respectively, the stable and unstable parts of the energies of relevant phases. The black circle 
denotes the ground-state energy of the disordered phase (DP). 

Fig. 5. Phase boundaries in the A — T c plane for the doubly decorated square Ising-Hciscnbcrg lattice 
(q = 4) when the exchange interaction Ji is changed. FP, QAP and DP denote the ferromagnetic 
(or ferrimagnetic), quantum antiferromagnctic and disordered phase, respectively. 

Fig. 6. Temperature dependences of the reduced specific heat of the doubly decorated square lattice (q = 4) 
for Ji/J = 0.5 and different values of the anisotropy parameter A. The dashed (full) lines represent 
the cases corresponding to the FP (QAP) ground-state phases, respectively. The insert shows the 
detail of the behavior when the anisotropy parameter A takes the value close to the boundary value 

A c = y/2. 

Fig. 7. Temperature variations of the correlation function qffc between Ising and Heisenberg atoms for 
q = 4, Ji/J = 0.5 and different values of the anisotropy parameter A. 

Fig. 8. The same as in Fig. 7 but for the correlation q^ z h between nearest-neighboring Heisenberg atoms. 

Fig. 9. The same as in Fig. 7 but for the correlation qf^ — qf y h between nearest-neighboring Heisenberg 
atoms. 
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